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IMPROVEMENTS IN RADIO POSITIONING SYSTEMS_ 

The present invention relates generally to improvements in radio positioning systems 
and their methods of operation and, more particularly, to methods and apparatus for 
5 simplifying the acquisition of data required in such systems. 

EP-A-0 303 371, the contents of which are hereby incorporated by reference, describes 
a radio navigation and tracking system which makes use of independent radio 
transmitters set up for other purposes. The signals from each transmitter, taken 
10 individually, are received by two receiving stations, one at a fixed and known location, 
and the other mounted on the mobile object whose position is to be determined. A 
representation of the signals received at one receiving station is sent via a link to a 
D processor at the other receiving station, where the received signals are compared to 

?J find their phase differences or time delays. Three such measurements, made on three 

y 15 widely spaced independent transmitters, are sufficient to determine the position of the 

Jp mobile receiver in two dimensions, i.e. its position on the ground. The phase or time 

offset between the master oscillators in the two receivers is also determined. 

j? "CURSOR", as the system described in EP-A-0 303 371 is known, is a radio 

! " 20 positioning system which can use the signals radiated by existing non-synchronised 

radio transmitters to locate the position of a portable receiver. Unlike some other 
systems which use the temporal coherence properties of networks of purpose-built 
synchronised transmitters, CURSOR makes use of the spatial coherence of the signals 
transmitted by single transmitters. In a further development (see EP-A-0 880 712 & 
25 WO-A-99/21028), the technology has been applied to find the position of a mobile 
phone handset in a GSM or other digital telephone system, and these are examples of 
an 'Enhanced Observed Time Difference' (E-OTD) method using the down-link 
signals radiated by the network of Base Transceiver Stations (BTS) of the telephone 
system. 
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In the digital mobile telephone application described in EP-A-0 880 712, the contents 
of which are hereby incorporated by reference, the signals from each BTS within 
range of the handset are received both by the handset itself and by a fixed nearby 
receiver, the Location Measurement Unit (LMU), whose position is accurately known. 
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Representations of the received signals are passed to a Mobile Location Centre (MLC) 
where they are compared in order to find the time difference between them. Figure 1 
shows the geometry of a standard two-dimensional system. The origin of Cartesian 
co-ordinates x and y is centred on the LMU positioned at O. The orientation of the 
5 axes is immaterial, but may conveniently be set so that the y axis lies along the north- 
south local map grid. The handset, R, is at vector position r with respect to the LMU 
position O. A BTS, A, is shown at vector position a. 

Consider first the signals from BTS A. The time difference, At a , measured between 
10 the signals received at R and O is given by 

Ar a =(|r-a|-|a|)/u+^, 

where v is the speed of the radio waves, £ is the clock time offset between the clocks 
15 in the receivers at R and O, and the vertical bars each side of vector quantities denote 
that it is the magnitude of the vectors which are used in the equation. The value of s 
represents the synchronisation error between the measurements made by the two 
receivers. Similarly, may be written for two other BTSs (B and C) at vector positions 
b and c (not shown): 

20 

Ar b =(|r-b|-|b|)/y+f , 

At c =(|r-c| - \c\)/u + s . ' (1) 

25 The values of At a , A/ b , At c , are measured by the methods disclosed in EP-A-0 880 712 
and the values of a, b, c, and u are known. Hence the equations (1) can be solved to 
find the position of the handset, r 5 together with the value of £ . 

In WO-A-99/21028, the contents of which are hereby incorporated by reference, it is 
30 described how these same time offsets can be measured using locally-created 
templates in a GSM telephone system as follows. Suppose that the handset R has 
recorded a short burst of the GSM signals from BTS A. Contained within that 
recording is the framing structure, synchronisation bursts and other 'given' data (or 
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predetermined values) which are a constant feature of those transmissions. The 
processor within the handset can create a matching template, based on the known 
structure of the network signals. Received signals can then be matched by the locally- 
generated template. When the template finds a match, the correlation peak at the 
position of best match corresponds to the time offset between the received signals and 
the local clock inside the handset. For the signals radiated by BTS A this measured 
time offset, Ar a i, is given by 

A/al =(l r " a IV^ + £ \ 3 

where a* is the time offset of the BTS transmissions and S\ is the time offset of the 
handset's internal clock, both relative to an imaginary universal 'absolute 1 clock. The 
signals from BTSs B and C may also be measured in the same way, giving 



15 A/bi =(|r - b\)/u + a h + s u 
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and 



At cl =(|r - c\)/v + Oc+ €\ . (2) 
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ljl The same measurements can also be made by the LMU, giving 



A/a2 =(|a|)/t> + a*+ £l , 
25 A/ b2 =(|b|)/u + «b+ ei , 

and 
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Ar c2 =(|c|)/t> + «c+ £ 2 , (3) 

where £ 2 is the time offset of the LMU's internal clock relative to the same imaginary 
universal absolute clock. Subtracting equations 3 from equations 2 gives 
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Af a = Af a i - A^a2= (|r - a| - \z\)/u + £ , 



A/ b = Afci - At b2 = (|r - b| - |b|)/u + e , 



and 



At c = At c \ - At c2 = (|r - c| - \c\)/u + ^ , 



(4) 
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where € = S\ - £ 2 . It will be noted that equations 4 are just like equations 1, and can be 
solved in the same way to find the position of the handset, r, and the value of s . 

It will be apparent that the CURSOR method as described above, in common with all 
other methods which use the signals from non-synchronised transmitters, requires a 
network of LMUs to be set up within the coverage area of the telephone system. 
These units act as reference points at which the unsynchronised signals radiated by the 
BTSs are measured for comparison with the same signals received by a handset. Each 
position measurement requires a match to be made between the signals received by the 
handset from a number of nearby BTSs, and signals received by an LMU from the 
same set of BTSs. In practice, it is often difficult to find a match using just one LMU, 
especially if the LMU network is sparse, since the handset may receive signals from 
BTSs not received by the LMU, and vice-versa. It is therefore necessary to combine 
the measurements from two or more LMUs. However, each new LMU brought into 
the calculation adds a further unknown clock time offset (s 2 , £3 etc.), each of which 
therefore requires an additional BTS measurement to provide the extra equation 
needed to solve for all the unknown quantities. 

One solution to this problem is presented in WO-A-99/21028 where it is shown how 
the LMU network can be synchronised. Referring to Figure 2, suppose that an 
adjacent pair of LMUs, Uj and U 2 , can see a common BTS. The positions of the 
LMUs and the BTS are all known, so a single measurement of the BTS signals by 
each LMU is sufficient to determine the clock time offset between the LMUs. For 
example, suppose that the distance from Ui to the BTS is s\ 9 and the distance from U2 
to the BTS is s 2 . Ui measures time offset A*i and U2 measures At 2 , given by 
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A/i = s\/v + a + S2\ , 

At 2 = s 2 /u + a + €22, (5) 

5 

where a is the time offset of the BTS transmissions, and s 2 \ and £22 are the time 
offsets of the LMU internal clocks in Ui and U 2 respectively. Subtracting the second 
equation from the first yields 

10 £21 - € 2 2 = Afi -Ar 2 + s\/u-s 2 /u, (6) 

which is the relative time offset of the clock in Ui with respect to that in U 2 . This 

~2 process may be repeated for a second pair of LMUs, say U 2 and U3, and another BTS 

03 whose signals can be received by both members of this second pair of LMUs. In this 

p 15 way a synchronisation map may be calculated, which provides the clock offsets of all 

% the LMU internal clocks relative to one of them adopted as a master 4 LMU network 

SI clock time'. Having established the LMU synchronisation map in this fashion, a 

O CURSOR position measurement can then include any number of LMUs without the 

T. penalty of adding an extra unknown time offset for every LMU, since the relative 

ru 

St 20 LMU time offsets are known. 

O 
H* 

The receivers discussed in the preceding paragraphs make measurements of time 
offsets. More generally, receivers can measure time offsets, phase offsets (which can 

be converted into time offsets with a modulo 360° ambiguity), frequency offsets or 
25 rates of change of frequency offsets. Though these measurements are of different 
quantities, the present invention is applied usefully to each of them as, when combined 
with similar measurements made by a second receiver, they can independently provide 
positional information. Positioning systems making use of these measurements are 
discussed in a related patent application (our reference MJB06427WO) filed 
30 simultaneously herewith. 
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The present invention teaches how the same advantages of an (effectively) 
synchronised LMU network may be obtained by setting up one or more 'virtual 
LMUs' in the network which act as interface nodes for the real LMUs. 

5 According to a first aspect of the invention, there is provided a method of generating a 
list of offsets in time, phase, frequency, or derivatives thereof, or their equivalents 
expressed as offsets in distance or derivatives thereof, of a plurality of transmission 
source signals, received at a given location, relative to a common reference, the 
method comprising 

10 (a) acquiring data from one or more receivers, the positions of which may 

be known or determined, the data from a receiver comprising offsets in 
time, phase, frequency, or derivatives thereof respectively of signals 
received from the transmission sources relative to a reference source in 
each receiver or to each other; and 

15 (b) combining the acquired data and calculating the list of offsets relative 

to the common reference. 

In practice the offsets from the list can be used in place of offsets obtained directly 
from the receiver or receivers. 

20 

The relative offsets in time, phase, frequency, or derivatives thereof, with respect to 
each other or with respect to a reference source, of the signals received by a first 
receiver from a plurality of the transmission sources may be represented by 
corresponding offsets or differences in the distances between the transmission sources 
25 and the first or second receivers. 

The invention also includes apparatus using the method above, the apparatus 
comprising 

(a) means for acquiring data from one or more receivers, the positions of 
30 which may be known or determined, the data from a receiver 

comprising offsets in time, phase, frequency, or derivatives thereof 
respectively of signals received from the transmission sources relative 
to a reference source in each receiver or to each other; and 
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(b) means for combining the acquired data and calculating the list of 
offsets relative to the common reference. 

In a method using techniques similar to or as described in EP-A-0 880 712, in place of 
5 offsets in time, phase, frequency, or derivatives thereof, or their equivalents expressed 
as offsets in distance or derivatives thereof, data representative of the received signals 
may be used, from which the offsets of signals received from the transmission sources 
relative to the reference source may be determined. 

10 Therefore, the invention also includes a method of generating a list of offsets in time, 
phase, frequency, or derivatives thereof, or their equivalents expressed as offsets in 
distance or derivatives thereof, of a plurality of transmission source signals, received 
O at a given location, relative to a common reference, the method comprising 

m (a) acquiring data from one or more receivers, the positions of which may 

JSf 15 be known or determined, the data from a receiver being representative 

*F ■ ■ of the received signals; 

Sj (b) determining from the acquired data the offsets in time, phase, 

p frequency, or derivatives thereof respectively of signals received from 

Hp the transmission sources relative to a reference source or to each other; 

rU 

%i 20 and 

~ (c) combining the offsets so determined and calculating the list of offsets 

relative to the common reference. 

The invention also includes apparatus for carrying out the method described 
25 immediately above, the apparatus comprising 

(a) means for acquiring data from one or more receivers, the positions of 
which may be known or determined, the data from a receiver being 
representative of the received signals; 

(b) means for determining from the acquired data the offsets in time, phase, 
30 frequency, or derivatives thereof respectively of signals received from 

the transmission sources relative to a reference source or to each other; 
and 

(c) means for combining the offsets so determined and calculating the list 
of offsets relative to the common reference. 
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A radio positioning method and system including either of the methods and apparatus 
defined above also form part of the present invention. 

The invention also includes apparatus (a 'virtual LMU') for carrying out either or both 
of these methods. The apparatus may include a computer (located anywhere 
convenient) and programmed to carry out the required process. Although the 
following description of a particular application of the invention concerns signals in a 
digital telephone network, it will be apparent that the invention is by no means 
restricted to this application but may be applied to any network of one or more 
transmitters, synchronised or unsynchronised, set up for any purpose. 

A virtual LMU includes a computer process which may run on any computer platform 
able to obtain data from real LMUs. Accordingly, a further aspect of the invention 
includes a method of calculating and maintaining a list of offsets in time, phase, 
frequency, or derivatives thereof, or their equivalents expressed as offsets in distance 
or derivatives thereof, of a plurality of transmission source signals, received at a given 
location, relative to a common reference. 

It is assumed that the network of BTSs is unsynchronised in that the transmission time 
offsets of the BTS signals bear no constant or known relationship to each other, but 
that nevertheless the BTS oscillators are quite stable, so that their instantaneous 
frequencies change only slowly with time. In these circumstances, it is possible to 
predict the currently received offset in time, phase, frequency, or derivatives thereof of 
the signals from a given BTS by a given real LMU from sufficiently-recent historical 
data. The real LMUs in the network make measurements of all the BTSs they can 
detect in a cyclic fashion, repeating the cycle every few seconds. They maintain these 
measurements in a stack, replacing the oldest measurements with the most recent. A 
linear or low-order polynomial fit to the measurements therefore provides a predictor 
for extrapolation into the near future, or for interpolation in the recent past. Assume 
that the BTS oscillators are sufficiently stable that reliable predictions can be made 
over a period of, say, ten minutes. Then, every few minutes, the virtual LMU 
(VLMU) contacts each real LMU and receives its predictors for the received offsets of 
the signals from all the BTSs in its measurement set. It is likely that many of the 
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BTSs will have been measured by more than one LMU, so the VLMU analyses the 
complete data set from all the real LMUs using well-known methods to determine 
both (a) the best values of the real LMU internal clock offsets in time, phase, 
frequency, or derivatives thereof, and hence (b) the received offsets in time, phase, 
5 frequency, or derivatives thereof of the signals from all the BTSs which would have 
been measured by a real LMU located at the assumed position of the VLMU and 
capable of receiving signals from every BTS. 

In the above description of the function of the VLMU, it should be understood that any 
10 or all of the LMUs could be replaced by other receivers, not necessarily fixed or at 
known positions, which have not been set up specifically as LMUs. For example, the 
data from a number of handsets could be used to determine frequency offsets if the 
handsets were stationary. Furthermore, it is shown in a related patent application (our 
reference MJB06427WO) filed simultaneously herewith, how the positions and 
15 velocities of handsets can be determined without the need for any LMUs at all. 

Particular advantages of using a VLMU in a network include the following: 
(a) a full match can be guaranteed between a handset's measurements and a single 
(virtual) LMU; 

20 (b) the VLMU process minimises the timing errors in individual LMU measurements; 

(c) the VLMU list is immediately available to the position calculation process, 
increasing the speed of computation; 

(d) when combined with the ideas described in a related patent application (our 
reference MJB06427WO) filed simultaneously herewith, a list of receive time 

25 offsets can be created in a network of BTSs where there are few if any real LMUs; 

(e) the VLMU provides, in effect, a synchronisation map of the real LMU network 
which network may then be used to monitor the BTS network and, in particular, to 
determine locations of newly-installed BTSs. 

30 One example of a method and apparatus according to the present invention will now 
be described with reference to the accompanying drawings, in which: 

Figure 1 shows the geometry of a CURSOR system as described in EP-A-0 880 712; 
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Figure 2 shows adjacent LMUs making measurements of a common BTS; 
Figure 3 shows a network of real and virtual LMUs in a system of the invention; 



5 Figure 4 illustrates a similar, simplified network; 



Figure 5 shows the positions of LMU sites and BTS sites in a real network; 
Figure 6 shows the same network with the addition of a Virtual LMU; 

10 

Figure 7 shows the same network with the real LMUs replaced by the Virtual LMU; 

Figure 8 illustrates, by way of a flow chart, the processing which takes place inside an 
LMU; 

15 

Figure 9 illustrates, by way of a flow chart, the processing which takes place inside the 
VLMU; 

Figure 1 0 shows a list of timing offsets generated in a real system such as that shown 
20 in the example of Figures 5 to 9; and 

Figure 1 1 shows a table of timing errors associated with LMUs and a virtual LMU in 
this example. 

25 By way of example, and with reference to Figure 3, the function of a virtual LMU is 
now described which determines timing offsets. 

Consider a network of TV real LMUs and M BTSs which includes a virtual LMU 
(VLMU). The position of the n^ LMU, U n , is represented by the vector u n and the 
30 position of the m 01 BTS, B m , is represented by the vector b m , both vectors being with 
respect to the same origin. Signals radiated by BTS m will be received by LMU n 
after a time lag, and the measurement of this time lag, A£ n m> is given by 
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Afnm = |u n -b m |/u + + a m ± cr r 



(7) 



where ^ is the clock time offset of LMU n, is the transmission time offset of BTS 
m, both with respect to an imaginary universal 'absolute' clock, and cr nm is an estimate 
of the error in the measurement of A* nm . The assumed position of the virtual LMU, V, 
is represented by the vector v. If the VLMU were able to receive the signals directly 
from BTS m without error then it would measure a receive time offset, J3 m , relative to 
the imaginary universal absolute clock, given by 



Substituting for Om in equation (7) using the value deduced from equation (8) gives 



Over the entire network of TV LMUs, all of the M BTSs are visible. Each individual 
LMU, however, will only see a few of them, but as long as there is significant overlap 
of visibility, it is possible to take the set of all At values and solve for values of Sn and 
p m . Hence the VLMU can calculate timings for any BTS as if the network of LMUs 
were synchronised, or as if just one LMU (the VLMU, V) was all that was needed to 
cover the entire BTS network. 

To illustrate this further, a simplified problem is shown and solved below using N = 2 
and M = 4, i.e. a network consisting of just 2 LMUs monitoring 4 BTSs (see Figure 
4). For simplicity it is chosen that s x = 0. This is allowed as the 'absolute' clock time 
is completely arbitrary and may, for example, be measured by the internal clock of 
LMU number one. (Note, however, that this choice does introduce an asymmetry into 
the solution in that the error associated with BTSs three and four are not the same.) 
The first LMU (Ui) can receive signals from BTS numbers one, two and three, but 
cannot receive the signals from BTS number four. The second LMU (U2) can receive 
the signals from BTSs one, two and four, but cannot see BTS number three. The 
equations can be written in matrix form as 



Pm = |v-b m |/u + CCm. 



(8) 



A/nm = |u n -b m |/u + £n + Pm~ |v-b m |/o ± CW 



(9) 
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or equivalently as A.x = b + Z, where Z is an unknown vector of the actual errors on 
each measurement. 

The standard technique known as 'least squares' postulates that the estimate for x 
which minimises Z is given by 



x = (A T WA)"' A T Wb, 



(11) 



where the symbol A indicates the transpose of the matrix A, and the matrix W is 
defined by 



w = 



'11 
0 

0 

0 

0 

0 



0 

1 



'12 

0 
0 
0 
0 



0 

1 



0 
0 
0 



0 
0 
0 

1 



'21 

0 



0 
0 
0 
0 

1 



'22 

0 



0 
0 
0 
0 
0 

1 



(12) 
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This particular example can be solved explicitly. For simplicity, it is assumed that all 
the values of cr nm are the same, and equal to cr. This gives the result 



fix = ( 3D, , + D 12 + D 2 \ - D 2 i )/4 
fi 2 = (3D ] 2 + D u +D 2 2-D 2l )/4 



±0.87<r, 
±0.87<7, 
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fa = ( 2D 24 + D u +D n - D 2X - D 22 )/2 
s 2 =(D 2l +D 22 -D u -D l2 )/2 



(13) 



± 1.41CT, 



± l.OOcr, 



where 



D nm = At nm - 

uni nnx 



K-bJ-lv-bJ 



(14) 
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Note that even in this simple case, when a BTS is seen by both LMUs, the errors in the 
calculated offsets are less than those in each of the measurements themselves. This is 
an important advantage of the virtual LMU method. 

The LMUs may also contain other means of synchronisation. For example, each real 
LMU could be connected to a GPS or other timing reference receiver, G, which serves 
to provide the common timing reference. In this case, the LMU network may be 
considered as synchronised already to this common timing reference (say GPS 
standard time), and then the VLMU need not solve for the individual values of e as 
these are already known. One advantage of using other means of synchronisation is 
that there is no longer a requirement for overlap of visibility of BTSs between adjacent 
LMUs. If each BTS site also carried an LMU, then that LMU would need only to be 
able to receive the (very strong) signals from its co-sited BTS transmitters), thereby 
simplifying the installation of the LMU antenna. 

The VLMU mode of operation described above may be called the 'pull mode' as it 
requires the VLMU to instigate data transfer to itself from every real LMU. It is also 
possible to have each real LMU continuously check the difference between its own 
prediction of the receive time offset from every BTS using the values calculated from 
the set of predictors last sent to the VLMU and the actual measured values. When any 
one of these differences exceeds a given value, the LMU can send its new predictor set 
to the VLMU. This mode of operation may be called the 'push mode'. The particular 
mode appropriate for a real system depends, amongst other things, on the stability of 
the BTS network. 
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An example of a prototype system constructed in accordance with the invention will 
now be described, in which timing measurements made by a number of LMUs are 
combined to create a list of receive time offsets for all the GSM BTSs in and around 
Cambridge, UK, as if they were observed by a single "Virtual" LMU. 

The positions of nine LMUs, Ui - U9 (shown as filled circles) and twenty three 
Cambridge area BTSs, B\ - B23 (shown by unfilled squares) are plotted in Figure 5 on 
an Ordinance Survey (OS) grid. Each LMU comprises (a) hardware including an 
internal clock, a GSM radio, a computer, and a telephone connection, and (b) software 
including a program to compile a list of receive time offsets. 

Figure 8 illustrates, by way of a flow chart, the main elements of the LMU software 
residing in each LMU, Ui - U9. Every few seconds the "scan cycle" program is 
entered at step Al. The program proceeds, in steps A2, A3, A6 and A7, to tune the 
GSM radio to each GSM transmission channel in turn and scan for BCCH signals. In 
the event that a BCCH signal is detected in step A3, the program calculates the receive 
time offset of the signal relative to its internal clock (in step A4) and updates its list of 
time offsets accordingly in step A5. The BCCH is also decoded to produce the ID of 
the BTS B n from which the signal is received. Figure 5 also shows lines between 
LMUs and the BTSs detected by them as a result of the scan cycle routine. 

Figure 6 shows the positions of the nine LMUs, the twenty three BTSs monitored by 
those LMUs and a Virtual LMU V (plotted as an unfilled circle). The Virtual LMU V 
comprises (a) hardware including a telephone connection to each real LMU and a 
computer, and (b) software including a program to compile a list of "virtual" receive 
time offsets. 

Figure 9 illustrates, by way of a flow chart, the main elements of the VLMU software. 
Every 4 minutes the "VLMU update" program is entered at step VI. The program 
proceeds, in steps V2, V3, V5 and V6, to connect to each of the real LMUs via the 
telephone connections. If the connection to a given LMU Ui - U9 succeeds, the 
program, in step V4, extracts that LMU's list of receive time offsets. On completion 
of the connection cycle, the program (in step V7) combines the data to produce a list 
of receive time offsets for all the BTSs monitored by the LMU network. An example 
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of this list is shown in Figure 10 where part of a table of receive time offsets generated 
by the VLMU is reproduced. The first column shows the BTS identifier and the 
figures in the table are in units of 1.85 microseconds. Figure 6 also shows lines 
between the VLMU V and each real LMU Ui - U9 from which the VLMU software is 
able to extract a list of time offsets. 

Once the "VLMU update" has been completed, the combined list of timings generated 
is equivalent to that which would have been observed had there been a single (real) 
LMU at the VLMU's location making timing measurements of every BTS in the 
network. Figure 7 illustrates this equivalence by showing the virtual monitoring of 
each BTS's timing measurements by the Virtual LMU (shown as lines between the 
VLMU V and the BTSs B, - B 23 ). 

Each timing measurement has an associated error which, in most cases, is smaller than 
the errors in the timing measurements made by the individual real LMUs. This is 
illustrated in the table of Figure 11, which shows part of the list produced by the 
VLMU during operation. The first column shows the identifier of each BTS. The 
columns headed 1 to 9 each relate to a particular one of the nine real LMUs. The 
figures in the table are errors in the timing measurements made by the real LMUs of 
the signals from the corresponding BTSs. An empty cell indicates that the LMU is 
unable to receive a BTS's signal. The column headed VLMU shows the result of 
combining the measurements using the method outlined above. 

The timings in the table of Figure 1 1 are in units of 1.85 microseconds. Note that the 
VLMU timing errors are generally smaller than those estimated for real LMU timings, 
confirming an advantage of the VLMU method in reducing errors. 
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